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ABSTRACT 

We present a new grid of model photospheres for the SDSS-III/APOGEE survey of stellar popu- 
lations of the Galaxy, calculated using the ATLAS9 and MARCS codes. New opacity distribution 
functions were generated to calculate ATLAS9 model photospheres. MARCS models were calculated 
based on opacity sampling techniques. The metallicity ([M/H]) spans from —5 to 1.5 for ATLAS and 
—2.5 to 0.5 for MARCS models. There are three main differences with respect to previous ATLAS9 
model grids: a new corrected H2O linelist, a wide range of carbon ([C/MD and a element [a/M] 
variations, and solar reference abundances from I Asplund et "ail (|2005| ). The added range of varying 
carbon and a element abundances also extends the previously calculated MARCS model grids. Al- 
together 1980 chemical compositions were used for the ATLAS9 grid, and 175 for the MARCS grid. 
Over 808 thousand ATLAS9 models were computed spanning temperatures from 3500 K to 30000 K 
and log g from to 5, where larger temperatures only have high gravities. The MARCS models span 
from 3500 K to 5500 K, and log g from to 5. All model atmospheres are publically available online. 



1. INTRODUCTION 

The Apache Point Observatory Galactic Evolution Ex- 
periment (APOGEE; Allende Prieto et al. 2008) is a 
large-scale, near-infrared, high-resolution spectroscopic 
survey of Galactic stars, and it is one of the four ex- 
periments in the Sloan Digital Sky Survey-Ill (SDSS-III; 
Eisenstein et al. 2011; Gunn et al. 2006, Aihara et al. 
2011). APOGEE will obtain high S/N, R-22,500 spectra 
for 100,000 stars in the Milky Way Galaxy, for which ac- 
curate chemical abundances, radial velocities, and phys- 
ical parameters will be determined. APOGEE data will 
shed new light on the formation of the Milky Way, as well 
as its chemical and dynamical evolution. To achieve its 
science goals, APOGEE needs to determine abundances 
for about 15 elements to an accuracy of 0.1 dex. To at- 
tain this precision, a large model photosphere database 
with up-to-date solar abundances is required. We chose 
to build the majority of APOGEE's model photosphere 
database o n ATLAS9 and MARCS calculations. 

ATLAS (IKuruczl fl979h is widely used as a univer- 
sal LTE 1— D plane-parallel atmosphere modeling code, 
which is freely available from Robert Kurucz's web- 
sitcfl ATLAS9 ([Kuruczl fl99l handles the line opac- 
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ity with the opacity distribution functions (ODF), which 
greatly simplifies and r educes t he computation tim e 
(IStrom fc Kuruc3 [19661: IKuruczl [20051: ICastellil [20051) . 
ATLAS uses the mixing-length scheme for convective en- 
ergy transport. It consists in pretabulating the line opac- 
ity as function of temperature and gas pressure in a given 
number of wavelength intervals which cover the whole 
wavelength range from far ultraviolet to far infrared. For 
computational reasons, in each interval the line opacities 
are rearranged according to strength rather than wave- 
length. For each selected metallicity and microturbulent 
velocity, an opacity distribution function table has to be 
computed. While the computation of the ODFs is very 
time consuming, extensive grids of model atmospheres 
and spectrophotometric energy distributions can be com- 
puted in a short time once the required ODF tables are 
available. 

ATLAS12 HKurucdl2005t ICastellill2005h uses the opac- 
ity sampling method (OS) to calculate the opacity at 
30,000 points. The high resolution synthetic spectrum 
at a selected re solution can then be o btained by run- 
ning SYNTHE dKurucz fc Avrettl ll98lD. More recently, 
iLester fc Neilsonl (|2008t ) have developed SATLAS__ODF 
and SATLAS._OS, the spherical version of both ATLAS9 
and ATLAS12, respectively. No extensive grids of models 
have been published up to now, neither with ATLAS12, 
nor with any of the two versions of S ATLAS. 

Instead, extensive grids of ATLAS9 ODF model at- 
mospheres for several m etallicities were calculated by 
iCastelli fc Kuruczl (f2003h . These grids are based on so- 
lar (or scaled sola r) abundance s from lGrevesse fc Sauvall 
(1998). Recently, [Kirbyl d2011[) provi ded a new ATLAS9 
grid, but he used abundances from lAnders fc Grevessel 
(1989). The calculations presented in this paper 
are based on t he mo re recent solar composition from 
lAsplund et al.l ()2005[ ). This updated abundance ta- 
ble required new ODFs and Rosseland mean opacity 
calculations as well. Abundances from lAsplund et al.1 



2 



were chosen instead of those from newer stud- 
ies (jAsplund et ah! I2009D to match the composition of 
the MARCS models described below, and those available 
from the MARCS website. 

The MARCS model atmospheres ([ Gustafsson et al.1 
[1975 IPlezetaJ] 119^1 iGustafsson et al.1 120081 ) were de- 
veloped and have been evolving in close connection with 
applications primarily to spectroscopic analyses of a 
wide range of late type stars with different properties. 
The models are one-dimensional plane-parallel or spher- 
ical, and computed in LTE assuming the mixing-length 
sch eme for convective e nergy transport, as formulated 
by iHenvev et al.1 ([1965). For luminous stars (giants), 
where the geometric depth of the photosphere is a non- 
negligible fraction of the stellar radius, the effects of the 
radial dilution of the energy transport and the depth- 
varying gravitational field is taken into account. Initially, 
spectral line opacities were economically treated by the 
ODF approximation, but later the more flexible and re- 
alistic opacity sampling scheme has been adopted. In 
the OS scheme, line opacities are directly tabulated for 
a large number of wavelength points (10 5 ) as a function 
of temperature and pressure. 

The shift in the MARCS code from using ODFs to the 
OS scheme avoided the sometimes unrealistic assumption 
that the line opacities of certain relative strengths within 
each ODF wavelength interval overlap in wavelength irre- 
spective of depth in the stellar atmosphere. This assump- 
tion was found to lead to systematically erroneous mod- 
els, in particular when plolyato mic molecules add i mpor- 
tant opacities to surface layers (|Ekberg et al.l l7986) . The 
current version of the MARCS code used for the present 
project and for the more extensive MARCS model atmo- 
sphere data_bas^3_was presen ted and described in de- 
tail by Gustafs son et al.l (|2008f ) . The model atmospheres 
presented in this paper adds large variaty in [C/M] and 
[a/M] abundances to the already existing grids by cov- 
ering these abundances systematically from -1 to +1 for 
each metallicity. 

Our main purpose is to update the previous ATLAS9 
grid and publish new MARCS models to provide a large 
composition range to use in the APOGEE survey and 
future precise abundance analysis projects. These new 
ATLAS models were calculated with a corrected H 2 
linelist. The abundances used f or the MARCS model s 
presented in this paper are f rom iGrevesse et al.l (l2007f) . 
which are nearly identical to lAsphmd et al.l (|2005h : the 
only significant differe nce is an abundance of scandium 
0.12 dex higher than in lAsplund et al.l(|2005D . The range 
of stellar parameters (T e ff, log g and [M/H]) spanned by 
the models covers most stellar types found in the Milky 
Way. 

This paper is organized as follows. In Section 2 we de- 
scribe the parameter range of our ODFs and model at- 
mospheres and give details of the calculation method of 
ATLAS9 we implemented. Section 3 contains the param- 
eter range and calculation procedure for MARCS models. 
In S ection 4 we compare MA RCS and ATLAS9 models 
with iCastelli fc Kuruca (|2003l ) , and illustrate how differ- 
ent C and a contents affect the atmosphere. Section 5 
contains the conclusions. The grid of ODFs and model 
atmospheres will be periodically updated in the future 

10 http://marcs.astro.uu.se/ 



[M/H] vs. [C/M] or [a/M] grid 
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Figure 1. The [C/M] or [a/M] content as a function of 
[M/H] of ATLAS9 models (filled circles, Table 1) and MARCS 
models (open circles, Table 4). Both [C/M] and [a/M] 
changes independently from each other, and the small steps 
in metallicities give altogether 1980 different compositions for 
the ATLAS9 models and 175 compositions for the MARCS 
models. The number of acceptable models may vary for 
each composition, for details see the ATLAS-APOGEE website 
flhttp: / /www. iac.es /proyecto /ATLAS- APOGEE / [ l . For missing 
metal-rich compositions of AlLASy models see table 2. 

and available onlinfl. 

2. ATLAS9 MODEL ATMOSPHERES 
2.1. Parameters 

The metallicity ([M/H]) of the grid varies from —5 to 
1.5 to cover the full range of chemical compositions and 
scaled to solar abundancetQ. For each of these solar 
scaled compositions we also vary the [C/M] and [a/M] 
abundances from —1.5 to 1 (Figure 1). 

ODFs and Rosseland opacity files were calculated with 
microturbulent velocities v t = 0, 1, 2, 4, 8 km s _1 , while 
the model atmospheres were produced only with vt = 2 
km s _1 . The metallicity grids were the same for all ef- 
fective temperatures, and the range can be seen in Table 
1. Some metal-rich compositions with high C but low 
a content were not calculated due to excessive compu- 
tation time; these are listed in Table 2. The a elements 
considered when varying [a/M] were the following: O, 
Ne, Mg, Si, S, Ca, and Ti. The temperature and gravity 
parameter grid for each composition and spectral type is 
given in Table 3. The T e g — log g distribution is plotted 
in Figure 2. Extreme metal-poor and metal-rich compo- 
sitions were also included. 

All the ATLAS codes use atomic and molecular linelists 
made available by Kurucz on a series ofCDROM(|Kuruci 
1990). They can now be found at the Kurucz websiteQ 
Th e molecula r line l ists f or TiO and H2O were p r ovide d 
by iSchwenkd (fl998l ) and iPartridge fc Schwenkel (fl997h . 
respectively, and reformatted by Kurucz in ATLAS for- 
mat. These are also available for download at the Kurucz 
website. For these models, we used the same linelists as 

11 http://www.iac.es/proyecto/ATLAS-APOGEE/ 

12 [M/H] means any element with Z > 2 and [M/H] = 
log 10 (N M /N H )* - log 10 (N M /N H )®, where N Fc and N H are the 
number of the desired element and hydrogen nuclei per unit vol- 
ume, respectively 

13 http://kurucz.harvard.edu/LINELIST/ 
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Abundance Parameters of ATLAS9 models 
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Figure 2. The gravity (log g) as a function of effective temper- 
ature (T eff ) of ATLAS9 and MARCS models calculated for each 
composition. Acceptable models are denoted by filled circles, while 
models not acceptable (and missing) are denoted by open circles for 
the solar composition. Models with C/O > 1.7 and T eff < 4000 K 
are not published (see Section 4.2). 



ICastelli fc Kuruczl (|2003j ). except for H2O, for which a 
new Kurucz release of the IPartridge fc Schwenkd (|1997t ) 
data was adopted 11 !. 

The solar ref e rence abundance table was adopted from 
lAsplund et al.l ()2005l ). Convection was turned on with 
the mixing length parameter set to l/H p = 1.25, but the 
convective overshooting was turned off. All the models 
have the same 72 layers from log tr oss — —6.875 to 2, 
where the step is log TR nss = 0.125. These pa rameters 
remained the same as ICastelli fc KuTucil (l200l for easy- 
comparison. 

All computations were performed on the Diodo cluster 
at the Instituto Astrofisico de Canarias. Diodo consists 
of 1 master node and 19 compute nodes, for a total of 80 
cores and 256 GB of RAM, communicating through two 
independent Gigabit Ethernet networks. 16 of the com- 
pute nodes host 2 Intel Xeon 3.20 GHz EM64T processors 
each, with 4 GB of RAM (2GB per core); the remain- 
ing 3 compute nodes host each 16 Intel Xeon (E7340) 
2.40 GHz EM64T processors, with 64 GB of RAM (4GB 
per core). On this cluster, about 3 months of computer 
time was required for ODF and model atmosphere cal- 
culations, using all 80 processors. 

2.2. Calculation Method 

Two separate scripts were developed, one for the ODF 
and Rosscland opacity calculations, and one for the AT- 
LAS9 calculations. The ODF and Rosseland opacity 
calculations follow e d exa ctly t he procedure de scribed by 
ICastelli fc Kuruczl (|2003l ). and iCastellil §005) using the 
DFSYNTHE code for the ODF, KAPPA9 code for the 
Rosseland opacity, and ATLAS9 for the model atmop- 
shere calculations (|Sbordond 12004 [2005). These codes 
were compiled in Linux with the Intel Fortran compiler 
version 11.1. 

Our algorithm sets up the initial starti n g mod els from 
the grid provided by ICastelli fc Kuruczl (pOOl FI. The 
algorithm chooses the model that has the closest com- 

14 http://kurucz.harvard.edu/MOLECULES/H20/h2ofastfix.bin 

15 http:/ /wwwuser. oat. ts.astro.it/castelli/ 





Min Max 


Step 


M/H] 


-5 -3.5 


0.5 


C/M] 


-1 1 


0.5 


a/M] 


-1 1 


0.5 


M/HJ 


-3 0.5 


0.25 


C/M] 


-1.5 1 


0.25 


a/M] 


-1.5 1 


0.25 


M/HJ 


1 1.5 


0.5 


C/M] 


-1.5 1 


0.5 


a/M] 


-1.5 1 


0.5 




Table 2 






Missing compositions of ATLAS9 models 




[M/H] 


[C/M] 


[a/M] 


1 


1 


-1.5 


1 


1 


-1 


1.5 


0.5 


-1.5 


1.5 


1 


-1.5 


1.5 


1 


-1 


1.5 


1 


-0.5 


1.5 


1 






Table 3 

Model Atmosphere Parameters of ATLAS9 models 



Spectral Type 


T c ft 


Teff 


Teff 


log 9 


log g 


log g 




min 


max 


step 


min 


max 


step 


M, N, R, K, G 


3500 


6000 


250 





5 


0.5 


F 


6250 


8000 


250 


1 


5 


0.5 


A 


8250 


12000 


250 


2 


5 


0.5 


B 


12500 


20000 


500 


3 


5 


0.5 


B, O 


21000 


30000 


1000 


4 


5 


0.5 



position, effective temperature, and log g to the desired 
output, and an initial ATLAS9 model is calculated. The 
result must be checked to see whether the output model 
satisfies the convergence parameters provided by the user 
for each layer in the model atmosphere. These param- 
eters were set to 1% for the flux or 10% for the flux 
derivative errors after 30 iterations in each run, as rec- 
ommended in the ATLAS cookboolj 1 ^!. A model is con- 
sidered converged if the convergence parameters satisfy 
these criteria in all depths. 

We then determined that an atmospheric model is ac- 
ceptable if one of the following criteria is satisfied: 1. 
the model has converged through the whole atmosphere, 
2. no more than 1 non-converged layer exists between 
log tr oss = —4 and log tr oss = 1. The model is allowed 
to have other non-converged layers for log tr oss < —4. 
A model was considered unacceptable in all other cases. 
We used only log tr oss >= —4 to log tr oss = 1 to check 
the convergence, because most of the lines in the optical 
and H band form in this region. In case the output was 
not acceptable, we restarted the calculation using more 
iterations. In case of a run with unacceptable output, we 
selected a starting model that had a different log g from 

16 http://atmos.obspm.fr/index.php/documentation 
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Table 4 

Abundance Parameters of MARCS models 





Min 






Max 




Step 


M/H] 


-2.5 






0.5 




0.5 


C/M] 


-1 






1 




0.5 


a/M] 


-1 






1 




0.5 






Table 5 








Model Atmosphere 


Parameters of MARCS models 




Spectral Type 


T B ft 


Teff 


Toff 


log g 


log g 


log g 




min 


max 


step 


min 


max 


step 


M. N. R, K a 


3500 


4000 


100 





3 


0.5 


K, G a 


4250 


5500 


250 





3 


0.5 


M, N, R, K b 


3500 


4000 


100 


3.5 


5 


0.5 


K, G b 


4250 


5500 


250 


3.5 


5 


0.5 



a Spherical atmospheres 

b Plane-parallel atmospheres 



the initial starting model and used it to restart the cal- 
culation. Then the previously described convergence test 
was performed and more restarts were done if it was nec- 
essary. If the output remained unacceptable, the effective 
temperature of the output model was changed by 10, 50, 
and 100 K. If the output of any of these runs was accept- 
able, then it was used as an input model to calculate the 
atmosphere with the original effective temperature. 

To test our model atmospheres, we replicated the 
Rosse land opacity calculations by iCastelli &: Kuruc3 
(2003). For this we used Castelli's scripts without 
any modification to calculate the PDF , and Rosse- 
land opacities for the iGrevesse fc Sauval (|1998l) abun- 
dances. These calculations concluded in perfect agree- 
ment within numerical precision. We then attempted 
to reproduce model atmospheres with the same param- 
eters found on Castelli's website with our scripts us- 
ing t he ODFs and Rosseland opacity files generated 
with IGrevesse fc Sauval (|1998fl abundances. This test 
also concluded with near perfect agreement with only 
0.1— 0.2 K maximum differences coming from the differ- 
ent version of compilers used. 

3. MARCS MODEL ATMOSPHERES 
3.1. Parameters 

Models were computed for seven overall metallicities, 
[M/H] from —2.5 to 0.5, with a step size of 0.5 dex. For 
each of these seven overall [M/H] mixtures, 25 combina- 
tions of modified carbon and a element abundances were 
adopted: The modifications to the logarithmic C and a 
abundances are —1, —0.5, 0, 0.5, and 1 dex. This format 
resulted in a total of 175 subgrids with unique chemical 
compositions (Table 4). The a elements in MARCS are 
O, Ne, Mg, Si, S, Ar, Ca, and Ti. This composition 
scheme is exactly the same as the previous models found 
on the MARCS website. The systematic a abundance 
changes were chosen to overlap the scheme (Section 2.1) 
used in the ATLAS9 model calculations. 

For each of these abundance subgrids, models with 12 
values of effective temperature from 3500 to 5500 K and 
1 1 values of logarithmic surface gravities from to 5 were 



computed (see Figure 2 and Table 5). Models with log- 
arithmic surface gravities lower than 3.5 (giants) were 
computed in spherical geometry and with a microturbu- 
lence parameter of 2 km s , while the remaining (dwarf) 
models adopted v t = 1 km s _1 and plane-parallel geom- 
etry. In the end, 86% of the 23,140 models converged 
satisfactorily. Convergence was particularly poor for cool 
dwarfs that are simultaneously a rich and carbon poor. 

Details about the atomic and molec ular line lists used 
by th e MARCS code are given by [Gustafsso n et al.1 
(2008). In a number of instances, they are different from 
those used by the ATLAS c ode, as for H 2 (|Barber et al.1 
l20M ). and TiO (|PlezHl998f l. 

3.2. Calculation Method 

Similarly to the ATLAS9 ODF calculations, a new 
metallicity MARCS subgrid is started with the summa- 
tion of an OS table of atomic line opacities for the rel- 
evant abundance mixture and microturbulence parame- 
ter. Since line opacity data is needed for many atoms 
and first ions, it saves time to add the opacities of the 
individual species into one file before a set of models is 
computed. This file contains a table with over 10 5 wave- 
length points with line opacities relevant to the equation 
of state for 306 combinations of temperature and damp- 
ing pressure Pq . The damping pressure is used as a proxy 
for the pressure that broadens metal lines by collisions 
with neutral atoms in cool stars. Pq is the pressure of HI 
with the addition of the polarizability corrected pressures 
of neutral helium and H 2 (see Eq. 33 of Gustafsson et 
al. 2008). For each molecule, in contrast, the line opac- 
ities are given in one table for the same wavelength set, 
but as a function of temperature and microturbulence 
parameter. 

The method used for the MARCS mod e l calc ulations 
is described in detail in iGustafsson et all (|2008l ). They 
start with the generation of a simplified starting model 
assuming a grey opacity. Physical parameters and their 
derivatives are computed and the model structure is it- 
erated in a multidimensional Newton-Raphson scheme 
until the flux through each depth layer corresponds to 
the prescribed effective temperature. The models usu- 
ally converged after 4 to 8 iterations; convergence also 
requires that the maximum temperature correction in 
any depth point is below 1.5 K during two consecutive 
iterations. Occasionally convergence takes longer, and 
some models do not converge at all. A converged model 
with similar model parameters is then identified as a new 
starting model. This approach is often successful, but 
some models do not converge, which leaves vacancies in 
the model grid. 

4. DISCUSSION 

Changes in the chemical composition can have a num- 
ber of effects on the calculated atmospheres. These ef- 
fects are mainly related to cither changes in opacity or 
changes in the equation of state. The main effect of an 
increased line opacity in late-type stars is either cooling 
or warming (depending on the opacity and its wavelength 
dependence) of the ou ter layers and also bac k-warming 
of the innermost ones (|Gustafsson et al1 ll975). Changes 
in the equation of state are mainly variations in the mean 
molecular weight for abundant elements, changes in the 
number of free electrons for elements that are important 
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Figure 3. Examples of ATLAS9 and MARCS model atmospheres 
with T cff = 4000 K and log g = 1 (red line) and T cff = 5500 K 
and log g = 4 (blue line). The model from the APOGEE grid is 
denoted by a solid line, the corresponding ATLAS9 model from 
ICastelli fc Kuruczl H2003) is denoted by a dotted line, the MARCS 
model is denoted by dashed line. The panels show the mass column 
(m), temperature (T), gas pressure (P), electron number density 
(N e ) as a function of the logarithm of optical depth (log t), and 
the relative difference between the ATL AS9 (A K ), MARCS (Am ) 
calculations presented in this paper and Castclli & Kurucz ( 2003) 
(Ok). 

electron donors, or other more intricate changes related 
to chemical equilibrium through molecule formation. In 
the next two sections, we give examples of MARCS and 
ATLAS9 models from our grid, and illustrate briefly 
the changes in the model atmospheres related to large 
changes in C and a elements. 

4.1. Comparing our MARCS and ATLAS9 models to 
the Castelli- Kurucz grid 

Figure 3 illustrates examples of the atmospheric struc- 
tures for two models with solar composition, one for 
T off = 4000 K and log g = 1 (red line) and the 
other for T cS = 5500 K and log g = 4 (blue line). 
There are four double panels showing the parame- 
ter dependence with Rosseland optical depth of the 
mass column (top- left panel), temperature (top-right), 
gas pressure (bottom-left) and electron number density 
(bottom-right). The dashed, solid, and dotted lines 
correspond to the new APOGEE MARCS, APOGEE 
Kurucz, and earlier (NEWODF) Castelli-Kurucz mod- 



els l|Casteili fe Kuruczl [20031 ) . respectively. While the 
APOGEE Kurucz and MARCS m odels share essentially 
the same chemical composition (|Asplund et al.l 120051 : 
IGrevesse et al.ll2007j) , t he Castelli-Kurucz models use the 
solar mixture given by IGrevesse" fc Sauvall (fl998h . 

For both the cooler (red line) and the warmer model 
(blue line), we see good agreement between the MARCS 
and the ATLAS9 models presented in this paper. The 
new MARCS and ATLAS9 models show good agree- 
ment, the differences are modest, less than 1 percent, 
for the thermal structure, and less than 2—3 percent 
for the gas pressure and the electron density in the lay- 
ers where weak spectral lines and continuum form (r > 
0.01). These small differences are most likely related 
to the different equation of state implemented in the 
two codes, and can also be due to the fact that AT- 
LAS9 uses the ODF method, while MARCS uses the OS 
method. The differences between the two increase at 
temperatures lower than 4000 K due to the different H2O 
and TiO linelists used in the calculations. More impor- 
tantly, larger differences are present in the gas pressure 
and electron density in the ATLAS9 models compared 
to the Castelli-Kurucz ones for the cooler models. These 
differences must be related to the updated solar chemi- 
cal composition, as the Castelli-Kurucz models are also 
computed with a corrected H2O linelist. The most signif- 
icant change in the solar composition corresponds to the 
reduction in oxygen, nitrogen and carbon abundances, 
all of which decrease in the update. This causes a re- 
duction in the Rosseland opacity at a given temperature 
and gas pressure. The decreased Rosseland opacity leads 
to a subsequent increase in the total pressure, which is 
only partially compensated by a reduction in the electron 
d ensity. 

iGustafsson et al. (200(1) present some further compar- 
isons between MARCS m odels and ATLAS9 models of 
ICastelli fc Kunic3 ((20031) . Other comparisons between 
MARCS, ATLAS9, and PHOENIX mode l atmospheres 
and spectra were discussed by IPlezl (f20lTh . 

4.2. Changes in model structures related to chemical 
composition 

Figures 4 and 5 illustrate the changes in the model 
structures in response to changes in chemical composi- 
tion. The left-hand panels of Figure 4 (T eS = 4000 K 
model) and 5 (T c g = 5500 K) show the relative changes 
associated with variations of the carbon abundance from 
solar proportions. The ATLAS9 model for the warmer 
temperature has Vt = 2 km s" 1 , while the MARCS model 
has Vt = 1 km s _1 , but this difference does not affect the 
model structures significantly. The relative variations 
found are very similar for MARCS and ATLAS9, except 
in the outermost layers of the atmosphere. The changes 
in the thermal structure are modest and show a behavior 
symmetric to that found for oxygen, which strongly sug- 
gests that CO formation is the driver of the variation. 
CO is the most tightly bound molecule and consumes 
almost all of the free atoms of cither carbon or oxygen, 
whichever is less abundant, leaving the majority species 
(C I or O I) to form other molecules. If oxygen domi- 
nates, it produces a "normal star" (if cool, an M star), 
otherwise carbon is free to form many molecules with 
high opacities, making carbon stars with very different 
spectra. At T e g = 5500 K, CO formation is low, thus 



6 



0.05 



-0.05 



0.5 



CL 

< 



-0.5 



Model: T eff = 4000 K, log g = 1 , ATLAS9 vs. MARCS 




ATLAS [C/M] = -1 

MARCS [C/M] = -1 



I I I I I I I I I 



I I I I I I I I I I 




ATLAS [a/M] = -1 

1 — 

MARCS [a/M] = -1 

1 — 



I I I I I I I I I I 




I I I I I I 



JA 




-7 -6 -5 -4 -3 -2 -1 
lo 9 T Ross 



1 2-7 -6 -5 -4 -3 -2 -1 
lo 9 x Ross 



1 2 



0.05 



-0.05 



0.5 



CL 



-0.5 

-1.0 
1.5 

1.0 

0.5 



-0.5 
-1.0 



Model: T eff = 5500 K, log g = 4, ATLAS9 vs. MARCS 



ATLAS [C/M] = -1 

MARCS [C/M] = -1 



I I I I I I I I I 



I I I I I I I I I I 



ATLAS [a/M] = -1 

1 — 

MARCS [a/M] = -1 

1 — 



I I I I I I I I I I 



J I I I I I I I I I 




-6 -5 -4 -3 -2 -1 
lo 9 T Ross 



1 2-7 -6 -5 



-4 -3 -2 -1 
lo 9 x Ross 



1 2 



Figure 4. Examples of how different [C/M] and [a/M] content 
changes the temperature profile in the ATLAS9 and MARCS at- 
mospheres for T c ff = 4000 K, log g = 1. This figure shows the 
temperature (T), gas pressure (P), electron number density (N e ) 
and the relative difference (where 0.5 corresponds to 50%, 1 corre- 
sponds to 100%) of these parameters between [C/M] = [a/M] = 
and [C/M] = [a/M] = — 1, and 1 as a function of optical depth (log 
r). The symbols T , P , N e>0 correspond to [C/M] = [a/M] = 0. 

this molecule does not affect the chemical equilibrium 
contrary to what is seen in the cooler atmospheres. 

The right-hand panels show changes for the same 
model parameters and different a-element abundances. 
Large changes are visible for both models in the pres- 
sure and electron numbers. The significant differences 
in the pressure compared to the cooler model are due 
to the increased gravity. Increasing the abundance of a 
elements reduces pressure (despite the fact that electron 
pressure increases), because the electrons of the other 
main electron contributors change the continuum opac- 
ity. The a elements that contribute most to the total 
number of electrons in the ATLAS9 models are shown in 
Figure 6. For the T cS = 4000 K, log g = 1 model, the 
main electron contributors are Ca, Mg, Na, and Al in the 
outer layers and Mg, Si and Fe in the deeper layers. For 
the warmer models, this changes significantly, as more 
Fe and H is ionized, the overall number of electrons in- 
creases and the main electron contributors become Mg, 
Si, and Fe through most of the atmospheres. 

Internal tests showed that polyatomic carbon 
molecules (C2H2, C3) substantially change the structure 



Figure 5. Examples of how different [C/M] and [a/M] content 
changes the temperature profile in the ATLAS9 and MARCS at- 
mospheres for T c fj = 5500 K, log g = 4. For more details see the 
caption of Figure 6. 

of the atmosphere with C/O ratios higher 1.7, if they 
are included in the linelist. These molecules significantly 
inflate the atmosphere changing the temperature in the 
line-forming photospheric layers. Since these molecules 
were used neither in the ATLAS9, nor in the MARCS 
calculations, atmospheric structures with C/O > 1.7 and 
T c ff < 4000 K are not reliable and thus not published. 

5. CONCLUSIONS 

Most of the ATLAS9 models are fully converged above 
T ff = 5000 K. Convergence problems are visible only in 
the outermost layers for stars with T e g < 5000 K. The 
regions affected by convergence issues are limited to log 
tross < — 4. However, these un-converged layers on the 
top of the atmosphere at low temperatures do not affect 
most line profiles very significantly, because most lines 
form deeper in the atmosphere. A little over a million 
ATLAS9 atmospheres were calculated. 

Over 20,000 MARCS fully converged models with 
ranges in [C/M] and [a/M] have been produced. Models 
are available in spherical and plane-parallel cases. Con- 
vergence issues are also present for red dwarf stars, sim- 
ilarly to the ATLAS9 models. 

Examples were presented for both ATLAS9 and 
MARCS models for T eff = 4000 K and 5500 K, and com- 
pared to the Castelli-Kurucz grid. The new ATLAS9 and 



7 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 


0.6 
0.5 
0.4 
0.3 
0.2 
0.1 


0.6 
0.5 
0.4 
0.3 
0.2 
0.1 


0.6 
0.5 
0.4 
0.3 
0.2 
0.1 





= 4000 K, log g = 1 T eff = 5500 K, log g = 1 T eff = 5500 K, log g = 4 





I I I I I I I I 



I I I I I I I I 



I I I I I I I I 



j i i i i i i i_ 



I I I I I I I I 



Si (a element) 



I I I I I I I I 



Ca (a element) 



I I I I I I I I 



Fe — 
Na — 

Al 




j i i i i i i i_ 



I I I I I I I I 



I I I I I I I I 



I I I I I I I I 



j i i i i i i i_ 



-6 -5 -4 -3-2-1012 -6 -5 -4 -3-2-1012 -6 -5 -4 -3 -2 -1 

log t Ross log t Ross log t Ross 



1 2 



Figure 6. Examples of the largest electron contributors to 
the total electron numbers in the ATLAS9 atmospheres for 
[M/H] = [C/M]=[o/M] = 0. The figure shows the element elec- 
tron numbers relative to the total number of electrons for the two 
models used in previous figures. The T e g = 5500 K, log g = 1 
model is plotted in the middle section to show that the electron 
contributors do not change significantly compared to a model with 
the same temperature but higher gravity. 

MARCS models agree well for all temperatures, while the 
differences between these calculations and the Castelli- 
Kurucz grid arise from the updated abundance tables. 
We briefly illustrated the effects of decreased/increased 
carbon and a content on the structure of the atmo- 
spheres, which are very similar in the MARCS and AT- 
LAS9 models. The response to the carbon content is 
only different in the outermost layers due to increased 
CO in the atmosphere. The response to a elements is 
also almost the same in MARCS and ATLAS9 models. 
The increased a— element content has profound effects 
in electron numbers and pressure for both the giant and 
dwarf stars, which is related to the higher number of Mg, 
Si, and Ca ions. Carbon rich models with C/O > 1.7 
and T c ff < 4000 K are not published here, because poly- 
atomic carbon molecules not included in the linelists 
here significantly change the temperature structure in 
the photosphere. 

These model grids will be used as the primary database 
for the pipeline analysing the spectra from the APOGEE 
survey. The high-resolution model spectra used in the 
survey's atmospheric parameters and abundances deter- 



mination code will be built on the models presented in 
this paper. Both the ATLAS9 and MARCS models at- 
mospheres will be continuously updated with new com- 
positions as the APOGEE survey progresses. The cal- 
culated ODFs, Rosseland opacities available for v t = 0, 
1, 2, 4, 8 km s -1 , and the ATLAS9 model atmosphere 
files are available for v t = 2 km s _1 from the ATLAS- 
APOGEE websit«Q. The MARCS models are available 
with vt = 2 km s" 1 for the giant and v t = 1 km s _1 for 
the dwarf stars from the standard MARCS websit<{3- 
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